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Abstract 

We study non-Fermi liquid effects due to the exchange of unscreened magnetic gluons in the 

normal phase of high density QCD by using an effective field theory. A one-loop calculation gives 

the well known result that magnetic gluons lead to a logarithmic enhancement in the fermion self 

energy near the Fermi surface. The self energy is of the form S(o;) ~ ury log(cj), where oj is the 

■ energy of the fermion, 7 = 0(g 2 ), and g is the coupling constant. Using an analysis of the Dyson- 

. Schwinger equations we show that, in the weak coupling limit, this result is not modified by higher 

£f) \ order corrections even in the regime where the logarithm is large, 7 log (a;) ~ 1. We also show that 

in ' 

this result is consistent with the renormalization group equation in the high density effective field 

in : 

theory. 
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I. INTRODUCTION 



It is theoretically well established that dense quark matter is not a Fermi liquid. Attrac- 
tive interaction between pairs of quarks that are anti-symmetric in color cause an instability 
in the quark-quark scattering amplitude if the quark momenta lie on opposite sides of the 
Fermi surface This instability is resolved by the formation of a diquark con- 

densate which breaks the color gauge symmetry. It is also well known that this is not 
the only non-Fermi liquid effect in dense quark matter. Unscreened magnetic gluon ex- 
changes lead to a logarithmic singularity in the quark self-energy close to the Fermi surface 
flflt QflflflQ. This legate may ,ea d t o a breakdown of pe rturb a- 

tion theory in the normal phase of dense QCD at very low energies, u ~ fiexp(—c n fl/g 2 ), 
where \i is the chemical potential, g is the coupling constant, and c n fl = 9n 2 . 

This scale is exponentially small as compared to the scale of superconductivity, u> ~ 
^iexp(— Cbcs/g), where Cb cs = 3tt 2 /^/2 [12]. Understanding non- Fermi liquid effects in the 
normal phase of quark matter is nevertheless important. First of all, understanding the 
normal phase is necessary in order to put calculations in the superfluid phase on a solid 
footing. Also, in order to establish the possible existence of a superconducting phase of 
quark matter from the observation of neutron stars we have to compute the properties of 
both the normal and the superconducting phase. And finally, non-Fermi liquid effects may 
play a role if the dominant superconducting phase is suppressed by non- zero quark masses, 
lepton chemical potentials, or magnetic fields. 

If electric charge neutrality is taken into account, then a non- zero strange quark mass 
leads to approximately equal differences between the Fermi momenta of strange and up as 
well as up and down quarks This implies that if the strange quark mass exceeds a 

critical value, only pairing between quarks of the same flavor is possible. Pairing between 
equal flavors requires order parameters with non-zero spin, and the corresponding gap s 
are suppressed by roughly two orders of magnitude compared to the spin zero gap |l4j . 
The gap can be further suppressed by a non-zero temperature or magnetic field. Finally, 
flavor symmetry breaking may lead to the appearance of gapless fermion modes even in the 



superconducting phase 15|,ll6j]. Depending on whether there is magnetic screening in this 



phase the gapless modes will also lead to interesting non-Fermi liquid effects. 

In this work we study non-Fermi liquid effects in QCD using the high density effective 
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theory [17|, [l8|, [ijj, 20]. The paper is organized as follows. In Sect. |TT] we discuss power 
counting in the high density effective theory. In Sects. IIHI and IIVI we study the Dyson- 
Schwinger equation for the quark self energy in the normal phase of dense quark matter. 
In Sect. |V]we consider the renormalization group equation for the quark propagator. We 
discuss some of the implications of our results in Sect. IV H Non- Fermi liquid effects due 
to unscreened gauge boson exchanges were first discussed by Holstein, Norton and Pincus 



in the case of a cold electron gas 
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23, 



24j . Similar effects due to d yna mical gauge 
fields in systems of strongly correlated electrons were studied by Polchinski 25|, Nayak and 
Wilczek [26], and others. 



II. HIGH DENSITY EFFECTIVE THEORY 



At high baryon density the relevant degrees of freedom are particle and hole excitations 
which move with the Fermi velocity v. Since the momentum p ~ v \i is large, typical soft 
scatterings cannot change the momentum by very much and the velocity is a ppr oximately 
conserved. An effective field theory of particles and holes in QCD is given by 



ipproxiu 

HQ 



(1) 



where t> M = {l,v}- The field if) v describes particles and holes with momenta p = fiv + I, 
where I <C We will write I = Iq + L + l± with h = v(l ■ v) and l± — I — In. In order to 
take into account all low energy degrees of freedom we have to cover the Fermi surface with 
patches labeled by the local Fermi velocity. 

Higher order terms in the effective lagrangian are suppressed by inverse powers of the 
chemical potential. There are two types of higher order corrections, operators that only 
involve fields in a given patch, and operators with four or more fermion fields that connect 
fields in different patches j^J. In order to understand the importance of higher order cor- 
rections we have to develop a power counting scheme for the high density effective theory. 
We first discuss a "naive" attempt to count powers of the small scale I. In the naive power 
counting we assume that v ■ D scales as I, ip v scales as Z 3 / 2 , scales as I, and every loop 
integral scales as Z 4 . We also assume that D±,v ■ D ~ I, where = (1, —v). In this case 
it is easy to see that a general diagram with V4 vertices of scaling dimension k scales as I s 
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FIG. 1: Counting hard loops in the effective field theory. If all (soft) gluon lines are removed the 
remaining fermionic loops contain sums over the velocity index. 

with 

5 = A + J2V k (k-A). (2) 

k 

A general vertex is of the form 

ip a (v ■ D) b (v ■ D) c (D±) d (l/ n) e , (3) 

and has mass dimension 3a/2 + b + c+d — e = 4. Since k = 3a/2 + b + c + d and e > we have 
k — 4 > 0. This implies that the power counting is trivial: All diagrams constructed from 
the leading order lagrangian have the same scaling, all diagrams with higher order vertices 
are suppressed, and the degree of suppression is simply determined by the number and the 
scaling dimension of the vertices. 

Complications arise because not all loop diagrams scale as Z 4 . In fermion loops sums 
over patches and integrals over transverse momenta can combine to give integrals that are 
proportional to the surface area of the Fermi sphere, 

1 r d 2 l ± jj, 2 r dfl 

2^^J J2nf ~2^J 4^' [ ] 

These loop integrals scale as Z 2 , not Z 4 . In the following we will refer to loops that scale as Z 2 
as "hard loops" and loops that scale as Z 4 as "soft loops". In order to take this distinction 
into account we define and V^ 1 to be the number of soft and hard vertices of scaling 
dimension k. A vertex is called soft if it contains no fermion lines. In order to determine 
the Z counting of a general diagram in the effective theory we remove all gluon lines from 
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the graph, see Fig. ^ We denote the number of connected pieces of the remaining graph 
by Nc- Using Euler identities for both the initial and the reduced graph we find that the 
diagram scales as I s with 

5 = J2[(k~ i)V k S + (* - 2 - h)Vf[ + Eq + A — 2N C . (5) 

k 

Here, fk denotes the number of fermion fields in a hard vertex, and Eq is the number of 
external quark lines. We observe that in general the scaling dimension 5 still increases with 
the number of higher order vertices, but now there are two important exceptions. 

First we observe that the power counting for hard vertices is modified by a factor that 
counts the number of fermion lines in the vertex. It is easy to see that four-fermion operators 
without extra derivatives are leading order (k — 2 — = 0), but terms with more than four 



fermion fields, or extra derivatives, are suppressed. This result is 
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amiliar from the effective 
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field theory analysis of theories with short range interactions 

The second observation is that the number of fermion loops that become disconnected 
if soft gluons are removed, Nc, reduces the power 5. Each disconnected loop contains at 
least one power of the coupling constant, g, for every soft vertex. As a result, fermion 
loop insertions in gluon n-point functions spoil the power counting if the gluon momenta 
satisfy I ~ gfi. This implies that for I < g\i the high density effective theory becomes 



non-perturbative and fermion loops in gluon n-point functions have to be resummed. 
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'his 



resummation corresponds to the familiar hard dense loop (HDL) resummation 
Note, however, that in the high density effective theory we do not perform a hard dense 
loop resummation of Green functions with external fermion lines. 

There is a simple generating functional for hard dense loops in gluon n-point functions 
which is given by j^] 

r - m2 V r a — r b (a 
Lhdl - — —2^ -r ; D yz^rf), V°) 

where m 2 = Nfg 2 fi 2 / (47r 2 ) is the dynamical gluon mass and the sum over patches corre- 
sponds to an average over the direction of v. For momenta I < g\i we have to add the HDL 
generating functional to the HDET effective action, Chdet — > C-hdet + C-hdl- This means 
that we work with hard dense loop resummed gluon propagators and vertices. If we analyze 
the low energy behavior in the vicinity of a generic point on the Fermi surface then there is 
no double counting involved, because we no longer have to consider sums over patches. It is 
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a) b) 





FIG. 2: Fig. a) shows the dominant contribution to the fermion self energy in the high density 
effective theory. The solid square denotes an insertion of the gluon self energy, see Fig. b). The 
dot denotes the free quark-gluon vertex. 

interesting to note that the velocity index on the quark field acts like a flavor label, and that 
the diagrams selected by the large \x limit are the diagrams of the large Nf approximation. 

The hard dense loop action describes static screening of electric fields and dynamic screen- 
ing of magnetic modes. Since there is no screening of static magnetic fields low energy gluon 
exchanges are dominated by magnetic modes. The resummed transverse gauge boson prop- 
agator is given by 

&0 - k 2 + lT)\ko\/\k\ 

where r\ = |m 2 and we have assumed that \ko\ < \k\. We observe that the gluon propagator 
becomes large in the regime k ~ (rjko) 1 ^ 3 ^> k . This implies that the power counting for 
very low energy gluons has to be modified. Landau damped gluons satisfy the scaling laws 
k° ~ I and \k\ ~ Z 1//3 . As we shall see in the next section, this modified scaling relation has 
important consequences for the structure of the fermion self energy. 



III. DYSON-SCHWINGER EQUATION 

The leading logarithmic term in the one-loop fermion self energy in the high density 
effective theory is determined by the Feynman diagram shown in Fig. El We find 



£(-,0 = |^iog(£), (8) 



where u = Iq and A is a cutoff. This result implies that for u ~ A exp(— 9n 2 /g 2 ) the 
perturbative correction is of order 1, and higher order terms of the form g 2n \og n (A / u) 
may have to be included. In order to study this problem we consider the Dyson-Schwinger 
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equation for the fermion self energy 

- <S(p) = - / ^f a »S(p + k)T b v D%{k). (9) 

Here, — iE(p) = tS' _1 (p) — S^ 1 ^) is the fermion self energy, D ab u is the gluon propagator, 
and T° is the quark-gluon vertex function. Following the arguments given in the previous 
section we consider the contribution from transverse gauge bosons only and use the HDL 
resummed transverse propagator. We will use the free quark-gluon vertex r" = gVi\ a /2. We 
shall solve the Dyson-Schwinger equation for the quark propagator but we will not solve a 
self-consistency equation for the gluon propagator or the quark-gluon vertex. We will justify 
these assumptions below. 

We note that the infrared divergence in the fermion self energy depends only on the 
energy and not on the momentum of the quark. Fermion momenta scale as l p ~ I while 
gluon momenta scale as \k\ ~ Z 1 / 3 . As a consequence we can neglect the dependence on the 
external fermion momentum and we will assume that the quark self energy is a function of 
the energy only. The Dyson-Schwinger equation is 

•w \ 2^ f d A k l-(v-k) 2 1 

- «S(po) = g C F - — -r : H — , (10) 

y "J (2ir) 4 Po + k -l p+k + i;(p + k ) k 2 -k 2 + i V \k \/\k\ 

where l p = v • p — pi and C F = (iV c 2 — l)/(2N c ). After analytic continuation to euclidean 
space we have 

2 f dk 4 f k 2 dk f 1 . 1 — x 2 I 



S(p 4 ) = g 2 C F / / — — / dx 



2tt J (2tt) 2 J-i i{p 4 + k 4 ) -l p -kx + iX(p 4 + k 4 ) kl + tf + r]^ 

(11) 

— * —* 

where l k = v ■ k = kx and we have defined £(p 4 ) = T lE (p 4 ) = —iH(p ). Because l p <C \k\ we 
can neglect the dependence on l p . The angular integration can be carried out analytically, 
leaving 

s(pi) = ^Lf dk J kdk ( fa±^±%±M!±j! arctan ( k \ 



4tt 3 J J \ k 2 \p 4 + k 4 + E(p 4 + k 4 ) / 

Pa + h + S(p 4 + A: 4 ) \ 1 

fc J k 2 + k 2 +r] \M> [ > 

This expression will be analyzed numerically in the next section. In order to derive an 
approximate analytic expression we note that the gluon propagator becomes large in the 
regime k ~ (rjk 4 ) l l z ^> k 4 . Therefore, we can in first approximation neglect the k\ term in 



the gluon propagator. In addition to that, we can approximate 1 — x 2 ~ 1 in the numerator 
of equ. (JTTJ). We get 

2 f dk^ f kdk 



S(p 4 ) = 2C F g 2 / — / - — - arctan - — ^. (13) 

The non-analytic contribution to the self-energy can be extracted from 

k ] i + s'fa + h) , s 

(k 4 +p A + J:(h+p 4 )) 2 + k 2 j P + r/M ' 1 ] 

Only the first term in the curly brackets has a logarithmic singularity in the limit p^ — > 0. 

We get 

which is equal to the result of the one-loop calculation. This implies that as long as g is small 
the one-loop self energy is a solution of the Dyson- Schwinger equation even in the regime 
g 2 log(A/p 4 ) ~ 1. It also means that there are no contributions of the form g 2n [log(A/p 4 )] n 
with n > 1. A similar conclusion was reached by Polchinski in his analysis of the spinon 
gauge theory in 2+1 dimension j^ . 

We finally return to the question of including vertex corrections in the Dyson-Schwinger 
equation. Brown et al. showed that vertex corrections are not logarithmically enhanced 
except in a small kinematic window where the collinear momentum transfer, AL, is much 

nn 

smaller than the energy transfer, Au; [8j, 20]. This conclusion is unchanged if self energy 
corrections to the propagator are included. The regime Al p Au; does not contribute to 
the Dyson-Schwinger equation at leading order in the coupling constant. As a consequence 
vertex corrections do not have to be included at leading order in the weak coupling limit. 
This result is analogous to Migdal's theorem for the electron-phonon interaction j3jj. We 
should emphasize, however, that the physical picture that underlies Migdal's theorem is 
quite different from what happens in the QCD case. 

We have also checked that fermion self energy insertions do not modify the gluon self 
energy, Fig. |2b, at leading order in the coupling constant. This result is related to the fact 
that the leading term in the self energy depends only on the energy, u, and not on the 
momentum, l p , of the quark. 
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TABLE I: Characteristic scales in dense quark matter. We compare the dynamical screening mass 
m, the BCS scale uJbcsi an d the scale of non- Fermi liquid effects ui n fl. We also show the one-loop 
running coupling constant a s evaluated at the chemical potential fj,. 

IV. NUMERICAL ANALYSIS 

In the previous section we presented analytic arguments which suggest that the infrared 
enhancement in the fermion self energy is one- loop exact in the weak coupling limit. In 
this section we shall strengthen these arguments by performing a numerical study of the 
Dyson-Schwinger equation (|12|1. This will also provide an estimate of the size of higher 
order corrections at non-asymptotic densities. 

Let us first give numerical estimates for the relevant scales. We shall assume that the 
strong coupling constant a s is given by the Nf = 3 one-loop running coupling constant 
evaluated at the scale fi. In Table H] we compare the dynamical screening scale m, the scale 
of superconductivity Ub cs ~ bofig~ 5 exp(— Cb cs /g), and the scale of non- Fermi liquid effects 
Unfi ~ m exp(— c n fl/g 2 ). We use Cb cs = 3n 2 /\^2 and c n fl = 9n 2 given above as well as 
bo = 5127r 4 exp(— (tt 2 + 4)/8). Even if the chemical potential is very large the screening scale 
is close to the Fermi energy. The scale of superfluidity varies very little whereas the scale of 
non-Fermi liquid effects becomes extremely small if the chemical is large. 

We have solved the Dyson-Schwinger equation for two different values of the chemical 
potential, a low value of /i = 0.5 GeV relevant for the physics of neutron stars and an 
asymptotically large value of fi = 100 GeV. At the high scale perturbation theory is expected 
to be applicable but at the lower scale the coupling is not small and the usefulness of 
perturbation theory is in doubt. It was shown, however, that the numerical solution of the 
gap equation in the superconducting phase is close to the asymptotic solution even if the 
coupling is not small, and that the gap on the Fermi surface is quite consistent with values 
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1 



fj, = 100 GeV 




100 



P4 [GeV] 

FIG. 3: Numerical solution of the Dyson- Schwinger equation for the fermion self energy. We show 
the euclidean self energy £(234) as a function of p^. The two sets of curves correspond to two 
different chemical potentials, \i = 0.5 GeV and [i = 100 GeV. The solid lines show the numerical 
solution, the dashed curves show the analytic result for the leading logarithm, and the dotted 
curves show a power-like self energy £(j>4) ~ p^K/p^) 1 . 

obtained from phenomenological models [3^. We have chosen the cutoff of the effective 
theory to be equal to the dynamical screening scale m. 

Our results are shown in Fig. The solid line shows the numerical solution and the 
dashed line shows the analytic result for the leading logarithm, equ. (|8j). For comparison we 
also show a self energy function that scales as a fractional power of energy, £(p4) = p^A/p^) 1 
with 7 = g 2 /(9n 2 ). This behavior was proposed by Boyanovsky and de Vega on the basis of 
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FIG. 4: Comparison between the numerical solution of the Dyson- Schwinger equation and the 
analytic leading log result. The solid lines show the ratio of the numerical solution over the leading 
log term for two different chemical potentials. The dotted lines show the same ratio for the power- 
like self energy H(p 4 ) ~ p 4 (A/p 4 ) 7 '. 

a renormalization group study [^] . Since 

P( (f) , , Pi { 1+7 l. 6 (f) + o(flo 5 '(f))} (16) 

this functional form starts to deviate substantially from the one-loop result for energies 
below the non- Fermi liquid scale, see Table HI 

For very small values of the energy, p4 C A, we find excellent agreement between the 
numerical results and the leading log expression. This is seen even more clearly in Fig. 0] 
where we show the ratio of the numerical solution over the leading logarithm. Significant 
deviations only occur for energies near the cutoff, but some higher order corrections are 
still present for energies one or two orders of magnitude below the cutoff. Our results show 
no evidence for higher order terms of the form g 2n log n (p4) with n > 1 and support the 
arguments given in the previous section. We conclude that in the limit of weak coupling 
and low energy the self energy is given by equ. (JBJ) even in the regime where the logarithm 
is large. 
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V. RENORMALIZATION GROUP 



Boyanovsky and Vega argued, on the basis of a renormalization group analysis, that the 
inverse quark propagator is of the form S~ 1 (uj, I) = uj{uj / A) 7 — I where 7 = g 2 /(9ir 2 ) in the 
weak coupling limit. A similar result in the QED case was presented by Gan and Wong 
I23I ] . These arguments contradict the solution of the Dyson-Schwinger equation presented in 
the preceding two sections. In order to clarify the situation we consider the renormalization 
group equation for the two-point function in the high density effective theory. We shall 
consider the leading order lagrangian 

C = i>l(u- v F l) if) v + gvpi^v ■ A^ v + . . . , (17) 

where we have explicitly included the Fermi velocity Vp. This is necessary because at finite 
density Lorentz invariance is broken and we need separate wave function renormalization 
factors for the energy and momentum dependent terms in the action. The one-loop fermion 
self energy is given by 



S(«,0 = Vfwlog(-), (18) 



g 2 v F /A" 
-tttu log - 

where we have neglected terms that do not contain logarithms. Equ. ()18|) has a logarithmic 
divergence in the effective field theory. This divergence can be removed by adding a counter- 
term to the lagrangian. We define the bare lagrangian C = C + C ct 

£ = V4 ( ZuJ - ZZ F v F l) ip v + Z g gv F tplv ■ Atp v 

= ipl iV (u - vq.fI) ?po,v + goVo.Fipl^v ■ Aip 0>v , (19) 

as well as the bare fields and coupling constants 

V>o,« = Z 1/2 ip v , v 0:F = Zpvp, g = -777T9- (20) 

Zj Zip 

Equ. (|18J) implies that at one-loop order Z ~ log(A) and that ZZp = 1. We showed in the 
previous section that, in the kinematic regime of interest, there is no logarithmic divergence 
in the quark-gluon vertex. As a consequence we can take Z g — 1. Equ. (|T%|) also suggests, 
however, that the effective coupling constant is 

9 2 vf (oU 

a = ^T' (21) 
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which acquires an anomalous dimension because of the scaling of the Fermi velocity. The 
bare and renormalized Green functions are related by 

G^\u t ,v , F h,a ) = Z n / 2 G (n \u t ,v F h,a), (22) 

where n denotes the number of external fermion fields, (a>j, k) are the external energies and 
momenta. Differentiating this relation with respect to A gives the renormalization group 
equation 

|a^ + p{ a )lL - lF{a)l~ + | T («)} h, a) = 0, (23) 

where we have defined the beta function and the anomalous dimensions 

x .da . . t d\ogZ , . k d\ogZ F . . 

^ = dA' 7(a) = A ^~' lF{a)=A ^k^- (24) 

In deriving equ. (}2~3*j) we have used the fact that G^ n ' depends on the Fermi velocity only 
through a and Vpli. At one-loop order we have 

/3(a) = -y F (a)a, 7(a) = -y F (a) = — . (25) 

We note that the beta function is positive and the effective theory is infrared free. This 
means that the perturbative analysis of the low energy behavior is reliable. The fact that 
the effective coupling is weak at low energy is related to the fact that the Fermi velocity 
vanishes as the quasi-particle energy goes to zero. We have checked that the one- loop 
anomalous dimensions do not depend on the gauge parameter in a generalized Coulomb 
gauge. In general, of course, there is no reason to expect the anomalous dimensions to 
be gauge invariant. Only physical properties of the solutions of the renormalization group 
equation, such as the quasi-particle properties, are gauge invariant. 

Boyanovsky and de Vega solved the renormalization group equation under the assumption 
that the beta function vanishes |2|. In this case we have 



4- 



S(w, I, a) = 0, (26) 



where we have used 7^ = —7. We observe that the inverse propagator satisfies the renor- 
malization group equation {A<9/(<9A) + y[ld/(dl) — 1]}S ,_1 (c<j, /, a) = 0. It is easy to see that 
this equation is solved by 

S-\lu,1) =uo (-) -v F l, (27) 
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where we have imposed the boundary condition S^ 1 ^ = A, I) = uo — vpl. Equ. (|27j) is the 
result of Boyanovsky and Vega. We showed, however, that the beta function does not vanish. 
The complete renormalization group equation is 



a d_ + l d_ + 1 
da dl 



S(w, I, a) 



0. 



(28) 



The propagator with the one-loop self energy included 



5- 1 (^,/)=^(l+ 7 log 



(29) 



is a solution of the complete renormalization group equation. We can also study the possible 
presence of higher order terms of the form S^ 1 ~ a n log n (A/uj). Consider the ansatz 



S 1 (UJ,1) =J2 a k 



a 10 



log 



— Vpl. 



(30) 



Inserting this ansatz into the renormalization group equation we obtain a^ + i = a^i{k — 
l)/(k + l), where we have used the one-loop anomalous dimension 7(a) = 7i<x This shows 
that dk = for k > 2 and terms of order a 2 log 2 (A/co>) or higher are absent. 



VI. SUMMARY AND DISCUSSION 



We have studied non-Fermi liquid effects due to unscreened transverse gauge boson ex- 
changes in the normal phase of high density QCD. We find that if the coupling is weak 
the fermion self energy is given by 70; log(A/c<j) with 7 = g 2 /(9ir 2 ). This result is reliable 
even if 7log(A/a>) ~ 1. We established this result using two different methods, the Dyson- 
Schwinger equation and the renormalization group. In the context of the Dyson-Schwinger 
equation the absence of higher order corrections is a consequence of the special kinematics 
of ungapped fermions interacting with Landau damped gluons. In the kinematic regime of 
interest the right hand side of the Dyson-Schwinger equation is independent of the fermion 
self energy. As a consequence, there is no difference between the one-loop result and the 
self-consistent solution. In the context of the renormalization group the absence of higher 
order terms follows from the relations 7 = — 7^ and (3 = —7^0 [24]. The relation between 
the anomalous dimension of the fermion field and the Fermi velocity is again due to the 
special kinematics. The relation between the beta function and the anomalous dimension of 
the Fermi velocity is a consequence of gauge invariance. 
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The weak coupling result implies that the quark propagator has a cut rather than a pole, 
and that the naive quasi-particle description breaks down. The spectral density is given by 

P (uj) = ^ 5 . (31) 

Ki+7iog(A/^))-/r + 

For non-zero momentum / this is approximately a Breit-Wigner distribution, but the wave 
function normalization and Fermi velocity vanish as / — > 0. An important consequence of 
the breakdown of Fermi liquid theory is an anomalous term in the specific heat. Ipp et 
al. showed that jll| 

C~ = 7 Cf ee log (J) = N f {Nl - 1)^ log , (32) 

where C£ ree = N c Nffi 2 T/3. They also computed the argument of the logarithm as well 
as terms that include fractional powers T 5//3 and T 7 / 3 . Our results suggest that equ. (|3^|) 
is reliable even if g 2 log(A/T) ^> 1. The anomalous term in the fermion self energy does 
not lead to an anomalous term in the thermodynamic potential Q at T = 0. The two- loop 
contribution to Q is infrared finite. Instead, this graph has an ultraviolet divergence in the 



effective field theory. This means that t 
in the microscopic theory. The result is 



re thermodynamic potential has to be determined 



In the superconducting phase the infrared enhancement in the fermion self energy is cutoff for 
energies less than the gap. Because A ~ exp(— Cb cs /g) the correction to the self energy never 
exceeds 7log(A) ~ 0(g). As a consequence non-Fermi liquid effects in the normal phase 
do not qualitatively alter the superconducting phase of QCD, but they give a correction to 
the gap which is enhanced by one power of 1/g relative to its naive order in the coupling 
constant. This correction reduces the gap by a factor exp[— (it 2 + 4)(AT C — 1)/16] ~ 0.18, 



where we have set N r 
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34 



35 



36|. 



In this paper we did not study Green functions with more than two external fermion 
lines. It is not clear whether more complicated n-point functions exhibit additional infrared 
divergences. It would be interesting, for example, to study the propagation of zero sound 
in the normal phase of dense quark matter. With regard to the physics of neutron stars it 
would also be interesting to study the thermal conductivity as well as the neutrino emissivity 
and opacity. 
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